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The aim of the present paper is twofold. The first goal is to show that high resolution angle- 
resolved photoelectron spectra from cuprates indicate the presence of string-like excitations of the 
quasihole excitation in a quantum antiferromagnet. In order to compare with the experimental 
intensity plots, we calculate the spectral function of the t — J and the t — t' — t" — J models within 
the self-consistent Born approximation for widely accepted values of the parameters of these models. 
The main features of the high resolution photoelectron spectra are in general agreement with the 
results based on the above models and can be understood by considering not only the lowest energy 
quasiparticle peak but, also, the higher energy string-like excitations of the hole. These features are: 
(a) the energy-momentum dispersion of the string-like excitations, (b) the momentum dependence 
of the spectral weight of the quasiparticle peak, and (c) the way in which the spectral weight of the 
lowest energy quasiparticle peak is transfered to higher energy string excitations, including the fact 
that it vanishes near the F point. The second goal of the present paper is to make the case that a 
proper analysis of both the numerical results obtained from these models and of the experimental 
results, suggests a theoretical picture for the internal structure of the hole-quasiparticle and a string- 
exchange pairing mechanism due to strong antiferromagnetic correlations among the background 
spins. We find that, using a simple model in which the Hilbert space is restricted to states of 
only unbroken strings attached to the holes and in which the holes are connected with unbroken 
strings, we can provide a good quantitative description of the most accurate numerical results. We 
find that the holes experience an effective interaction due to a string-exchange mechanism which is 
characterized by a rather large string tension and this can provide pairing energy scales much larger 
that those suggested by spin-fluctuation mediated pairing models. In addition, it is argued that 
such string-exchange interaction tends to bind holes in a bound state with the d x 2_ y 2 symmetry. 

PACS numbers: 71.10.-w,71.10.Fd,71.27.+a,74.72.-h,79.60.-i 



I. INTRODUCTION 

Angle-resolved photoelectron spectroscopy (ARPES) 
studies of the cuprates and related materials have been 
very useful in revealing the structure of the low energy 
quasiparticle excitations^^. In this paper, we will fo- 
cus on the most recent high resolution ARPES studies^ 
where an "anomalous" dispersion was identified at rel- 
atively higher energy than the well-known low energy 
quasiparticle band in the under-doped regime. We would 
like to discuss the case of very light doping where there 
are concrete predictions for the hole spectral function in 
a quantum antiferromagnet 5,6 . In this limit the problem 
has been extensively studied and from that ground and 
the above mentioned experimental studies a novel physi- 
cal picture emerges which consists of some already known 
features and some new ideas to be discussed here. 

The motion of a single hole inside a quantum antifer- 
romagnet has been extensively studied using a number 
of analytical, semi-analytical and numerical techniques. 
One of the main and well-known conclusions of these 
studies is that the hole becomes a rather well-defined 
low-energy quasiparticle with the minimum of the quasi- 
particle band at (±7r/2, ±7r/2), which is an unusual po- 
sition in the Brillouin zone for a minimum to occur, 
and, thus, at sufficiently light doping, a hole pocket de- 
velops near (±7r/2, ±7r/2). In addition, quantum spin- 



fluctuations allow the hole to move coherently as a po- 
laron because of the fact that they can erase the damage 
created by the hole and, vice versa, the hole through 
its motion which displaces reversed spins can eliminate 
some spin-deviations created in advance by quantum fluc- 
tuations; thus, the hole motion can lower the total en- 
ergy because it repairs some damage by "choosing" to 
move through paths of already reversed spins. Angle- 
resolved photoelectron studies of the undoped insulator 
have demonstrated^ that the minimum of the hole band 
is at (±7r/2, ±7r/2) and the bandwidth is in very good 
agreement to that predicted from the t — J model 5 when 
the experimentally determined valued of J ~ 0A3eV and 
the empirically known value of t ~ OAeV are used. The 
dispersion along the direction (0, 0) to (n, tt) as calcu- 
lated using the t — J model 5 does not agree with that 
measured by ARPESi and it requires an extension of the 
t — J model to add nearest neighbor (t 1 ) and next near- 
est neighbor (t") hopping to reproduce the experimental 
findings along this directionS&i&ii^ii&iiii 5 .. 

What is less well-known about the hole motion in a 
quantum antiferromagnet is what we call the string ex- 
citations which are internal excitations of the hole po- 
laron. These excitations arise because the hole along its 
path creates a string of spins which are reversed rela- 
tive to the local antiferromagnetically correlated back- 
ground. These excitations have been discussed by 
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Shraiman and Siggi a 16 ' 17 i 18 and investigated further by 
several author o 5 ' 20 i 21 ' 22 . In particular in Ref. [H, where 
a detailed calculation of the hole spectral function was 
reported, it was found that these excitations appear as 
rather sharp peaks in the hole spectral function and they 
compete with the low energy well-defined quasiparticlc 
peak to gain spectral weight as a function of momentum. 

The aim of the present paper is twofold. One of the 
goals of the present paper is to discuss that the re- 
cently reported high resolution angle resolved photoelec- 
tron studies from cuprates" strongly indicate the exis- 
tence of such string excitations in these materials (See 
also, Ref. A possible connection between such string 
excitations and certain features of the spectral function 
was made in past experimental studies^iiS.. In the present 
paper, a detailed analysis of their possible role in the 
spectral function is provided along with a study of the 
consequences of their possible presence in doped quan- 
tum antiferromagnets. We provide a rather transparent 
physical picture, through which we can understand sev- 
eral features of the ARPES results in terms of the same 
unified framework. The analysis provided in the present 
paper, by recalculating the hole spectral function using 
the t — J model 2 ^ and the t — t' — t" — J model in order 
to properly compare the results with the experimental 
data, makes a strong case for the existence of such string 
excitations in the real materials. 

The second important goal of the present paper is 
to present a new framework to understand the possible 
role of local strong antiferromagnetic correlations in the 
doped cuprate materials in the pairing mechanism. First, 
we introduce a model based on a reduced Hilbert space 
which contains only unbroken strings of overturned spins 
attached to the holes. It is demonstrated that this re- 
duced Hilbert space accurately captures the features of 
the quasiparticle dispersion and the higher energy string 
excitations. In addition, it is illustrated that using a sub- 
space of states in which a pair of holes is connected with 
unbroken strings leads to a string-mediated interaction 
which might give significantly stronger pairing mecha- 
nism than models of pairing due to exchange of spin- 
fluctuations. 

The paper is structured as follows: In Sec.[H]we present 
the results of the hole spectral function as calculated 
using the t — J and the t — t' — t" — J models pay- 
ing particular attention to the contribution and role of 
string excitations. In Sec. IIIH the analysis of the recent 
ARPES experimental results in order to reveal the con- 
tribution of the string excitations in the real materials 
is made. Sec. IIVI presents the consequences of our find- 
ing to our theoretical picture for the hole-quasiparticle in 
doped cuprates. In addition, a simple model is presented 
based on states of unbroken strings connected to the holes 
which gives a very good account of the numerical results 
of the t — J model. In addition, a pairing mechanism 
through string exchange interaction is presented in the 
same section. 



II. CALCULATION OF THE SPECTRAL 
FUNCTION 

In this section of the paper we will use the same 
method used in Ref. [H to calculate the spectral function 
of a hole in a quantum antiferromagnet using the t — J 
and t — t' — t" — J models with the aim to make a direct 
comparison with the recent ARPES results. 

In order to start the discussion, let us consider a quan- 
tum antiferromagnet doped with holes as described by 
the well-known t — J model 

n = -t ]T (4 a c ja + 4 a c ja ) + j z sis? 

<ij>,<7 <ij> 

+ ^tJ2(4s;+s 7s p (l) 

<ij> 

which is restricted to operate in a single occupied sub- 
space. Here c ia creates a spin a electron at lattice site i. 
The operators s^' y ' z are the three components of a spin- 
1/2 operator and sf — sf±isf. The first term takes into 
account the hole-hopping and the second and third are 
the direct and exchange terms in the Heisenberg antifer- 
romagnet. Just the first two terms form a Hamiltonian 
known as the t — J z model. The t — J model corresponds 
to the spacial case when J z = J±, The t — t' — t" — J 
model is obtained from the above t — J model by adding 
a hopping term t' which allows the hole to move along 
the diagonal of the square lattice and a t" term which 
moves the hole by two sites along the x direction (and 
— x direction) or by two sites along the y direction (and 
— y direction). The t — t' — t" — J model has been used 
in Refs. IMlMlIillllHH to rather successfully fit 
the photoemission data of Wells et alX 

In Ref. [1 following Ref. [23| an effective Hamiltonian was 
used where the Heisenberg spin interaction part of the 
t — J model was treated in spin- wave approximation and 
the hole-hopping term was linearized in spin-deviation 
operators. Within the non-crossing approximation, the 
Dyson's equation for the single-hole Green's function was 
solved self-consistently. When the calculation was car- 
ried out on small-size lattices a very good agreement was 
found with the exact diagonalization results for the spec- 
tral function. However, the advantage of the method is 
that it can be carried out on sufficiently large size lattice 
where finite-size effects are entirely absent^ i 24 i 25 . Hence, 
we will use this method to calculate the hole spectral 
function for the t — J and t — t' — t" — J models with 
the aim to make a direct comparison with the ARPES 
results. 

In Fig. Q] the hole spectral function is shown for the 
t — J model along the (0,0) to (ir/2,ir/2) direction us- 
ing the broadly accepted parameter values of t — OAeV 
and J ft — 0.3. The Dyson's equation was solved as in 
Ref.[5| using Lorentzian broadening with a width r\ = O.K. 
We note that while the width of the first peak depends 
strongly on the value of 77, as it is a 5 function peak, the 
other two (labeled II and III) and the higher peaks (not 



FIG. 1: The spectral function of a single hole in a quantum 
antiferromagnet for J/t = 0.3 along the (0,0) to (tt/2,tt/2) 
direction as a function of energy ui. 



clearly seen in the graph due to the choice of the scale) 
remain unchanged when we decrease the value of r\. This 
point has been extensively studied in Ref. [H and we have 
been able to reproduce it here. The first important point 
to be discussed in this paper again is the higher energy 
peaks labeled II and III in the spectral function. These 
peaks were studied in Ref. [H and it was concluded that 
they correspond to the string-like excitations which can 
be understood in the simple t — J z model. Before we go 
on in this paper we wish to briefly discuss the origin of 
these peaks. 

Let us first imagine a hole subject to quantum mechan- 
ics moving in a classical antiferromagnetic template such 
as the one conceptualized within the t — J z model. As 
the hole tries to move, it creates a string of spin devia- 
tions over the antiferromagnetic background and, thus, 
it experiences a linearly rising potentia l 16 ' 17 . As a result, 
if we ignore the small probability amplitude for the hole 
to retrace its pathi£, the hole is bound to its "birth" site. 
In addition, there are excited states ^ 20 ' 21 of the hole in- 
side the almost linearly rising potential which in the limit 
t/ J >> 1 can be approximated by the Airy functions and 
their corresponding eigenenergies are given by 

E n /t = e n + a n (J z /t) 2 / 3 . (2) 

Now, let us turn on the Heisenberg antiferromagnetic 
exchange coupling J±. This term creates or eliminates 
nearest-neighbor pairs of spin-deviations and, thus, it 
makes it possible for the hole to become free of its orig- 
inal "birth" site and hnd its "umbilical-cord" attached 
to a new "birth" origin. As shown in Fig. [5] this can 
either happen by eliminating the pair of spin-deviation 
at the beginning of the string involving the "birth" loca- 
tion of the hole (Fig. \2I.c) or anywhere along the string 
(Fig. [2^d), thus, producing two strings, one attached 
to the hole and the other becomes part of the back- 



(c) (d) 



FIG. 2: (a) A hole, shown by an open circle, in an anti- 
ferromagnetically aligned background (the spins of the Neel 
background are not shown in this figure for clarity; only the 
spin-deviations are shown by a solid circle) with no spin- 
deviations, (b) The hole has moved along the shown path 
by displacing the spin and creating the spin deviations shown 
with the black dots and thus creating the string shown, (c) 
The pair-spin-flip term of the Heisenberg antiferromagnetic 
exchange can eliminate a pair of nearest-neighbor spin devi- 
ation at the beginning of the string and, thus, shortening the 
string or (d) in the middle of the string thereby cutting the 
string into a string with the hole at one end and just a string 
of spin-deviations. 



ground spin-deviations, i.e., part of the fluid of the spin- 
deviations which always exist in a quantum antiferro- 
magnet due to zero temperature quantum fluctuations. 
In addition, to these processes there can be processes 
where the creation of pairs of spin-deviations adjacent 
to the hole position precedes the hole motion, in which 
case the hole experiences the string as a line of down- 
hill potential. In the limit of t/J± >> 1, however, the 
spin-relaxation time is much longer that the characteris- 
tic time for hole-hopping. Thinking in the same spirit as 
the Born-Oppenheimer approximation for the electronic 
motion relative to the slowly moving ions, one can con- 
vince himself that the approximate notion of a linearly 
rising potential exists. The important difference from 
the simpler t — J z model is that these spin fluctuations, 
no matter how slow they maybe, play the leading role 
in delocalizing the "birth" site from where the string is 
attached. This process is the process which gives rise to 
the hole band with a minimum at (7r/2,7r/2). We would 
like to stress, however, that a more correct qualitative 
picture is one where the hole is held by a string from a 
mobile birth-site. We will come back to this discussion 
in Sec. IIVI were we will examine further the role of these 
processes. 

In Ref. HI, the peaks II and III in Fig. [TJ have been 
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attributed to the internal excitations of the string analo- 
gous to the excitations of a particle in a linear potential. 
The energy where the peaks are located in the spectral 
function have been studied as a function of J/t (in the 
t — J model where Jj_ = J z = J) and it was found that 
they fit to the form given by Eq. [2] with J z replaced by 
J in a wide range of J/t. In addition, the width of these 
peaks was found to scale with J/t. In Sec. IIV[ the emer- 
gence of these peaks as higher energy excitations of the 
hole restricted in a Hilbcrt space of states made of only 
unbroken strings is discussed. 

In the present paper, one of the important aspects 
of Fig. [1] that we notice is the vanishing of the spec- 
tral weight of the low energy quasiparticle peak as we 
move away from the minimum of the band at (7r/2, 7r/2) 
toward the T (momentum (0,0)) point. In addition, no- 
tice that the spectral weight is transfered to the higher 
energy peaks which correspond to the string resonance 
excitations. 



I 




FIG. 3: The spectral function of a single hole in a quantum 
antiferromagnet for the t — t' — t" — J model with the same 
parameter values used in Ref. [l5l along the (0, 0) to (7r/2,7r/2) 
direction as a function of energy uj. 



In Fig. [3] the calculated hole spectral function for the 
t — t' — t"—J model with the same parameter values used 
in Ref. 15 (t = 0.35eV, t' = -O.UeV, t" = 0.08eV and 
J = O.UeV) is shown along the (0,0) to (n/2,w/2) cut. 
Again, in this case also, the Dyson's equation has been 
solved as in Ref. using Lorentzian broadening with a 
width ?y = O.lt. This figure should be compared with 
Fig. Q] obtained for the pure t — J model. All important 
aspects of the spectral function obtained with the simpler 
t — J model remain the same. The only difference is 
that as we approach the T point more spectral weight is 
transfered to the peak II instead of III. 

In Fig. [4] and Fig. [5l the calculated hole spectral func- 
tion along the (0, 0) to (it, 0) cut is shown for the t — J 
model and the t — t' — t" — J model respectively. The 
low energy quasiparticle excitation along this direction 




FIG. 4: The spectral function of a single hole in a quantum 
antiferromagnet for the t — J model along the (0, 0) to (n, 0) 
direction as a function of energy u>. 



attains it minimum at (7T, 0) in the t— J model (Fig. gj) 
where the peak acquires most of the spectral weight. This 
is not the case in the experiment. As a result the t — J 
model was extended to the t—t' — t" — J model precisely 
for this reason. Notice in Fig. [5] that the band minimum 
along this direction is no longer at (n, 0) but rather at 
(tt, /2, 0) where the spectral weight in maximum. Again, 
notice that the spectral weight as we approach the T 
point is transfered to the second peak due to the internal 
excitation of the string. 




FIG. 5: The spectral function of a single hole in a quantum 
antiferromagnet for the t — t' — t" — J model with the param- 
eters used in Ref.Ql (t = 0.35eF, t' = -0.12eV, t" = 0.08eV 
and J = O.lAeV) along the (0, 0) to (n, 0) direction as a func- 
tion of energy u. 
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III. COMPARISON WITH EXPERIMENT 




Energy (eV) Energy (eV) Energy (eV) 



FIG. 6: The experimentally determined spectral function with 
high resolution ARPES 3 shown on the left (part a) is com- 
pared with the theoretical spectral functions obtained from 
the t — J model* for J/t = 0.3 and t = OAeV along the path 
from (0,0) — > (7r/2,7r/2) shown in the middle (part b), and 
that obtained from the t — t' — t" — J model (part c) along 
the previously mentioned Brillouin Zone path. The solf-gray 
lines labeled 1,11 and III are guides to the eye. The reference 
energy has been shifted by a constant so that the energy of 
the quasiparticle peak at (tt/2, tt/2) to be the same as the 
experiment. 

The lowest energy peak at (tt/2, tt/2) in the theoretical 
calculation is very sharp while in the experimental data 
it carries a width of the order of a fraction of an eV. It 
is not clear what is the exact mechanism which causes 
the damping of the quasiparticle and there are reasons 
to believe that it may be due to the coupling to phonon 
excitations^,. Therefore, in order to compare with the 
experimental results we need to broaden the results of 
our calculation. We use a Gaussian broadening function, 
namely A^i/k, oj) — J A(k, (jj')G(lu — uj')dcu' , where G(lu — 
lu') = exp(-(uj-uj') 2 /e 2 ), using e = y^a = 0.125eV. We 
have obtained very similar results when we use a value 
of 7] = Q.125eV in the initial propagator namely in the 
iterative procedure to solve the Dyson's equation. The 
same amount of broadening has been used in Ref. [15| to 
compare the ARPES peak to the quasiparticle peak given 
by the t — t' — t" — J model. 

In Fig.[S]a comparison is presented between the exper- 
imental spectral function along the (0, 0) to (tt, tt) cut, as 
reported in Ref. H by means of high resolution ARPES 
(left or part (a)), and those obtained from the t— J model 
(part (b)) and the t — t' — t" — J model (part (c)). We 
would like to notice that both in the theoretical and the 
experimental results, as we move away from (7r/2,7r/2) 
toward (0, 0), the spectral weight of the low energy well- 
defined quasiparticle peak vanishes and it is gradually 
transfcrcd to higher energy excitations. Notice that in 
the comparison of Fig. [6] and in the comparison provided 
in Fig. [JJ and Fig. [51 if the energy scale of the theoreti- 
cal results (i.e., the value of t) is increased by about 25% 
(such that the theoretical peak at the T point at ~ 1.5eV 



is pushed to ~ 2eV) the agreement between theory and 
experiment becomes better, especially at higher energies 
near 2eV. However, at such high energy scales we expect 
other effects to become important and the description, in 
terms of a Zhang-Rice^ singlet within the simple t — J 
model, is also expected to fail. 

In Fig. [JJ we have produced an intensity plot along 
the (0,0) to (tt, tt) cut, in order to compare directly our 
results with those of Ref. U For the purpose of compari- 
son, we have added the experimental results of Ref. y as 
reported in their Fig. 2 (top of Fig.[jj. In the second in- 
tensity plot from the top, we present that obtained from 
the t — J model with the generally accepted parameter 
values mentioned earlier. The bottom intensity plot is 
that obtained from t — t' — t" — J model using the pa- 
rameter values used in Ref. Il5l The intense peaks at 
(it/2, tt/2) and near (0, 0) are present in both theoretical 
and experimental intensity plots. As we approach mo- 
mentum (0,0) the spectral weight is transfered to higher 
energy "string" states (II and III) as also was discussed 
in SecHIl This gradual transfer manifests itself as a more 
luminous path connecting the bright peaks at (7t/2, 7r/2) 
and (0,0). 

In a very recent ARPES study^ an anomalous disper- 
sion was identified and a second energy scale at around 
0.%eV was reported. In our calculation this energy scale 
is the average energy of the string excitation peaks II and 
III (measured from the lowest energy state at (tt/2, tt/2)). 
This becomes evident by comparing our Fig.[7Jwith Fig. 2 
of Ref. i. In Fig. E we present a theoretical color-coded 
intensity plot (left) along the (tt/2, tt/2) to (0,0) cut, as 
obtained from the t — J model to be compared with the 
three ARPES intensity plots (right part of Fig. [5J for un- 
derdoped, optimally doped and overdoped Bi2212 taken 
from in Fig. 1 of Ref. [4j. A very similar intensity plot 
was obtained for the case of the t — t' — t" — J model. 
First of all in making a comparison we must realize that 
the lowest energy quasiparticle states, which correspond 
to the lowest energy string states in the neighborhood 
of (tt/2, tt/2), are occupied by holes which form a Fermi 
sea. Namely, these occupied states will be entirely dark 
in the intensity plot and the states with highest intensity 
would be right on the intersection of the Fermi surface 
with the most luminous states of the theoretical inten- 
sity plot. Therefore, by assuming these adjustments the 
color-coded intensity plot of Fig. [5J excluding the bright 
spots around (tt/2, tt/2) (because at sufficient amount of 
doping the states around (7t/2, 7t/2) should be inside the 
Fermi sea and the bright spots should move just out- 
side the Fermi surface), qualitatively agrees with the ex- 
perimental intensity plots reported in Fig. 1 of Ref. 
These figures and the previous discussion clearly indi- 
cate that the process of spectral weight transfer to the 
higher string states is masked by low intensity and/or 
broadening and manifests itself in the intensity plot as an 
"anomalous" dispersion with the center of the anomaly 
close to (7r/4,7r/4) as in Fig. 1 of Ref. 0. We find it 
very interesting that these features persist all the way up 
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Momentum 



FIG. 7: The experimentally observed intensity plot reported 
in Fig. 2 of Ref. |3| (top) is compared with that obtained from 
the t — J model (middle) and from the t — t—t'—J model 
(bottom). The reference energy has been shifted by a constant 
so that the energy of the quasiparticle peak at (n/2,n/2) to 
be the same as the experiment. 



to the overdoped regime. Our conclusion, therefore, is 
that the spin correlations should be strongly antiferro- 
magnetic even in the overdoped cuprates and this point 
is discussed in the next section. 



IV. THE STRING MEDIATED INTERACTION 

A. Role of string states in the single-hole 
dispersion 

As we discussed the origin of the spectral function 
peaks is due to the string excitations and they provide 
a qualitative explanation of the vanishing of the spec- 
tral weight from the low energy quasiparticle peak by 
means of spectral weight shift to higher energy string ex- 




n/4 %I2 n/4 7t/2 n/4 nil it/4 ji/2 

Momentum * x (1/a) 



FIG. 8: Left: The theoretical color coded intensity plot along 
the cut indicated in the top-right inset, as obtained from the 
t — J model using the widely acceptable values for the parame- 
ters J/t = 0.3 and t = OAeV. Right: Three ARPES intensity 
plots for underdoped, optimally doped and overdoped Bi2212 
taken from Fig. 1 of Ref. Q for comparison. 



citations. In this part of the paper, motivated by these 
finding we will consider a new theoretical approach to un- 
derstand doped quantum antiferromagnets, one in which 
such string excitations play a fundamental role. 

First we notice that as the hole moves it does not sim- 
ply create spin deviations, it creates (or destroys) a spa- 
tially connected string of spin-deviations. Therefore, as 
a basis to span the Hilbert space instead of disconnected 
spin-deviations we will consider an over-complete basis 
which is spanned by all possible strings which define con- 
nected paths of spin-deviations as shown in Fig. [21 The 
spin fluctuation part of our Hamiltonian conserves spin 
and creates or annihilates pairs of spin deviations which 
can be regarded as strings of two sites. The only case in 
which a single spin-deviation is needed is when we have 
a state where a spin-deviation is next to the site where 
the hole is located. This state can be regarded as a hole 
attached to a string of length unity. 

The fundamental reason for the choice of this basis is 
that in the limit where t >> J, the quantum antiferro- 
magnet with just a small hole concentration should look 
like a fluid of strings of overturn spins created by the mo- 
tion of the holes with a lot of nearest-neighbor pairs of 
spin deviations (due to the spin conserving antiferromag- 
netic exchange) which can also be considered as strings 
of size two. In this regime, the average string size in the 
ground state should be larger than two. In addition, this 
basis describes better the short wavelength physics and 
it might provide a better description in the doped regime 
when the long-range spin correlations are destroyed, but 
the relatively short-range antiferromagnetic correlations 
are still present. 

To show the importance of string excitations, we diag- 
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(d) 



length 4 (solid curve). Notice that the agreement with 
the results obtained with the method discussed in Sec. [TT] 
is remarkably close. The absolute energy scale at k = is 
accurately reproduced and, in addition, notice the close 
agreement of the dispersion near (7r/2,7r/2). There is 
some disagreement near (jr, 0). However, this part of the 
dispersion is not described accurately by the full solution 
to the t — J model either and we needed to introduce the 
t' and t" terms which give the leading correction to the 
dispersion near this point. 

The previous calculation suggests the following varia- 
tional single-hole wavefunction, where the effect of back- 
ground spin-fluctuations above the Neel state is turned 



FIG. 9: Connected single-hole string states up to length 4. 
The hole is represented by the open circle and the filled circle 
stands for a spin deviation, (a) The hole in the Neel back- 
ground, (b) There are 4 states of string length unity, (c) 
There are 12 states of string length 2. (d) 36 states having 
string length 3 exist, (e) There are 108 states such the ones 
shown in Fig. (e) above with strings of length 4. 



where 



R 

n I 



(3) 



r R = i + EE (4) 



R+<5„ 



P * - P s> 



(5) 




10.0) 



FIG. 10: The hole dispersion as calculated with a Hilbert 
space made of only strings. The dispersion is shown along 
the Brillouin zone path (0, 0) -> (tt, 0) -» (tt, tt) -» (0, 0). The 
three dispersion relations are calculated with string states of 
up to a limiting length size which is indicated in the figure. 
The result is compared to those obtained from the full cal- 
culation of the dispersion curve from the spectral function 
discussed in Sec. [TT] 



onalized the t — J Hamiltonian in a restricted subspace of 
only states having a string connected to the hole. These 
states are shown in Fig. [JJ|a-e). They can be charac- 
terized by the length of the string and we have succes- 
sively included all such states with strings up to length 
4. In Fig. [TfJ we present the calculated dispersion rela- 
tions where the energy is obtained by including strings 
states of (a) up to string length 2 (dashed-dot solid line) 
(b) up to string length 3 (dashed curve) and (c) up to 



where Pr/,,5 = c R , a c R , , $ a is the hole-hopping oper- 
ator of the t — J Hamiltonian which creates the strings 
and ^^l({Si}) are variational parameters. The opera- 
tors Tj£({8i}) create strings of overturned spins of length 
m with the hole at R along the path R — > R + 5i — > 
R + 5*i + 5 2 -> ... -> R + A + 5 2 + ... + 5 m . The prime 
on the summation means that it is over paths such as 
<5, + i =/= —Si. The summation over m is over states of 
strings of length m (up to n) connected with the hole. 
The following state can be used for the state |0) in order 
to include the effect of spin-fluctuations 



10) = 



; ( -l)^>exp(-i 



i<j 



(6) 



where the sum is over all spin configurations and L(a) is 
the so-called Marshall phase 7 , namely L(a) is the num- 
ber of up spins in one sublattice in the configuration a. 
When Uij = the state |0) is the Neel state along the 
x-direction in the spin space. The variational wavefunc- 
tion given by Eq. [?] when Uij — is identical to that 
obtained from the diagonalization discussed previously 
(with results given in Fig. [TO")) . Furthermore the func- 
tion is the spin-spin correlation function and describes 
spin-fluctuations^ around the Neel ordered state. Vari- 
ational calculations for the Heisenberg antiferromagnet 
using the no-hole state given by Eq. [5] give very accurate 
results for the ground state energy, the staggered magne- 
tization and the excitations above the ground state via 
sum rules 7 . In Ref. l30l the above wavefunction with n = 2 
and the same as the one which describes the no-hole 
case accurately^ was used to carried a variational Monte 
Carlo calculation. The results of this calculation for the 
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single hole dispersion were accurate for J/t > 1. For 
J/t ~ 0.3, longer string states were required (n > 2), 
which is a conclusion very similar to the present findings 
presented in Fig. [TO] These longer string excitations were 
allowed in the Green's function Monte Carlo(GFMC) cal- 
culation of Ref. [H which iteratively projects the correct 
single hole state starting from such a wave function as 
an initial state. This calculation converges to the correct 
energy with only a few GFMC iterations which indicates 
that this wavefunction is a good starting point having 
relatively large overlap with the exact wavefunction. 



termining the energetics of the hole dynamics. One ex- 
pects that long wavelength spin fluctuations to be im- 
portant in the ordered quantum antiferromagnet. How- 
ever, when the antiferromagnetic order is destroyed as the 
doping level is increased, while the long wavelength spin- 
fluctuation excitations are not well-defined, the string ex- 
citations still remain and they can quantitatively account 
for the short wavelength dynamics. With these findings 
in mind we proceed to show that these string correlations 
give rise to a pairing mechanism which is characterized 
by a relatively strong pairing interaction, namely, sig- 
nificantly stronger that the effective interaction due to 
exchange of antiferromagnetic spin-fluctuations. 
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FIG. 11: The higher energy string excitations obtained by 
diagonalizing the t — J Hamiltonian in the subspace spanned 
by the 161 states illustrated in Fig. [9] using J/t = 0.3 and 
t — OAeV. The lowest energy state is labeled and it is 
the same as the dispersion presented in Fig. 1101 The lowest 
lying excited states are three nearly degenerate excited states 
labeled 1,2,3 (which form the group labeled (II)). Following 
a gap there is a group of states labeled (III) which includes 
the fourth state and the next six states, namely, 5-11 which 
are degenerate. In addition, it includes the 12th excited state 
which follows a mirror-like dispersion to the ground state. 

In Fig. QTJ we give the energy of the string excita- 
tions found by diagonalizing the t — J Hamiltonian in the 
subspace of the string states shown in Fig. [9] along the 
same Brillouin zone path, namely, from (0, 0) — > (n, 0) — > 
(-7T, 7r) — > (0, 0). These states form bands and we can iden- 
tify what we have labeled as II and III in the previous dis- 
cussion of Sec. [II] and Sec. IIIII Namely, the three nearly 
degenerate excited states labeled as 1,2 and 3 form the 
band II and the nine states labeled 4-12 form the group 
labeled III in our previous discussion of Sec. [Til Notice 
that their energy is the same to those of the peaks labeled 
II and III in Fig. [TJ and Fig.|H In addition, their disper- 
sion is very similar, namely, II is almost dispersionless 
and III along the (0,0) to (n/2,Tr/2) cut (Fig. [TJ moves 
closer to the peaks I and II near (0, 0) and further away 
from the peaks I and II near (7r/2, tt/2). 

Therefore, the conclusion of this study is that the 
short-range string correlations are very important in de- 



B. String- mediated pairing 

In this paper, inspired by the findings of the photoe- 
mission measurements and our suggestion that they pro- 
vide evidence for the existence of string excitations even 
at optimal doping, we suggest that the pairing mecha- 
nism found numerically in the t — J mode l 33 ' 34 is medi- 
ated by strings. 

The rather simple calculation of the previous subsec- 
tion suggests that states of unbroken strings are the main 
contributors to the single-hole dispersion. In Fig. [TJ] a 
pair of holes A and B is shown connected by a string of 
overturned spins in an otherwise idealized Neel ordered 
background with no spin deviation. First, in order to 
facilitate the discussion let us consider states with un- 
broken strings connecting the two holes. The two holes 
experience a potential linearly rising as a function of the 
length of the string between them. Namely, when either 
hole is moved in such a way that the distance between the 
holes is increased as shown in Fig. fl"27 b). it stretches the 
length of the string of overturned spins and the potential 
energy increases. On the other hand, each hole feels a 
down-hill potential toward each other because by mov- 
ing along the path, defined by their mutually connecting 
string, they repair the overturned spins as indicated in 
Fig. He). 

In addition, it can be demonstrated using Fig. [TJ] 
(a,b,d) that when there are such string correlations be- 
tween two holes, the entire string bipolaron can move 
through the lattice with an effective mass which is lighter 
than that of a single hole. In leading order the hole has to 
make two successive hops followed by a pair spin-flip in 
order for the hole to move through the lattice. This is of 
higher order than the process illustrated in Fig.[T27a,b.d) 
where first the hole A hops (transition from (a) to (b)) 
momentarily lengthening the string and this is followed 
by a hop of the second hole B (transition from (b) to (d)) 
which takes the down-hill path formed by the string. No- 
tice that the configuration (d) has the same length of the 
string as the starting configuration A, but the hole string 
bipolaron has moved and, thus, gains delocalization en- 
ergy. 

A variational wave function to describe such a pair of 
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FIG. 12: (a) A pair of holes (A) and (B) marked by o's are 
connected by a string of overturn spins, (b) For either hole to 
attempt to unbind itself from the other leads to lengthening 
the string which is energetically costly, (c) For either hole 
to attempt to hop closer to the other hole it is energetically 
favorable, (d) The string configuration of pair of holes (a) 
can be momentarily lengthened by allowing hole B to hop as 
shown in (b) and finally the length of the string is restored to 
original length size by allowing hole A to hop along the string 
as in configuration (d). 



holes of zero net momentum may be written as follows: 

ri ,r 2 

x [l+ J2 ^(^(ri ! r 2 ))y(r 1 ,r 2 )]f ri f r2 |0)(7) 

s(ri,r 2 ) 

where Y(ri, r 2 ) is a connected path connecting the sites 
ri and r 2 which are the location of the holes and 
<£(y(ri, r 2 ))) is a function of the path which defines the 
string and it is expected to be a decreasing function of 
the length of the path. The state |0) is the no- hole state 
given by Eq.[6]and the operators T ri 2 are given by Eq.[4] 
and they create the strings attached to each hole. The 
operator F(ri,r 2 ) is a product of hole-hop operators Py 
(defined after Eq. [5]) along the path connecting the two 
holes defined by F(ri,r 2 ). When this operator acts on 
the Neel state creates states, such as the state illustrated 
in Fig. [r2T a). where the holes located at ri and r 2 are 
connected with a string of overturned spins in an oth- 
erwise perfectly ordered Neel spin lattice. The function 
g(r 2 — rx) gives the symmetry of the two hole state and 
for d x 2_ y 2 symmetry it changes only its sign (keeping the 
same magnitude) for a n/2 rotation. 

A similar variational state built upon such a basis of 
string states which included the effect of background spin 
fluctuations was studied in Ref. and it was also used 
as the initial state for the Green's function Monte Carlo 



calculation. It was found that two holes bind in a state 
with the d x 2_ y 2 symmetry for J/t greater than about 
0.27. 

Let us restrict ourselves in the subspace where the two 
holes are connected by an unbroken string. This implies 
that the spin-fluctuation part of the t — J model gives 
zero contribution to matrix elements in this subspace. In 
order to illustrate that, let us consider a two- hole state 
with an unbroken string mediated between the two holes 
such as any state from Fig. [121 When applying the pair 
spin-flip operator of the t — J model on such a state, 
it will always produce a state with a broken string and, 
thus, it takes us outside the Hilbert space of states with 
all-connected strings. As a consequence the result of the 
calculation of the binding energy A of two holes in this re- 
stricted subspace, is very similar to that using the t — J z 
model with J z — J. Using a diagrammatic approach 
Chernyshev and Leung 36 found that there is strong bind- 
ing in p and d wave symmetry between two holes in the 
t — J z model. In addition, calculation of the quantity 
A has been carried out for the t — J z model using exact 
diagonalization 3 ^ on lattices with up to 50 sites and it 
was found that the two holes bind in a d x 2_ y 2 for J/t 
greater than ^0.2. 

We would like to go back and discuss the fact that the 
spin fluctuation operator can break the string connecting 
the two holes. These processes, however, compete with 
processes which lower the pair hole energy discussed pre- 
viously by allowing the pair to delocalize when the con- 
necting string is present and with other processes which 
accommodate the hole motion of each hole within the 
string. Namely, when one hole creates a string, another 
hole near the first hole, has increased probability to follow 
the string created by the first hole, thus, they are expe- 
riencing pairing correlations. Therefore, on variational 
grounds we believe that the pairs of holes, depending on 
the parameter regime, may have string-mediated pairing 
correlations. 

Here, we would like to discuss what happens when the 
level of doping increases. As it has been discussed earlier 
in this paper, in order for the string-like correlations to 
be effective, we do not need to require long range antifer- 
romagnetic(AF) order. This is because the cuprates are 
described by an intermediate parameter regime, where 
t/J ~ 3. In this intermediate regime, we found that 
strings of length only three lattice spacings can reproduce 
the single hole dispersion reasonably well (see Fig. [TU)) . 
Therefore, short-range string correlations only require 
that AF spin-spin correlations exist up to such rather 
short distance of two or three lattice spacings away from 
the hole. This requirement is believed to be true even at 
optimal doping. Hence, we believe that such short-range 
AF spin correlations, required for the string mediated in- 
teraction between pairs of holes, could still be present in 
the cuprates even at optimum doping. 
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V. CONCLUSIONS 

Our findings suggest a new framework to understand 
doped quantum antiferromagnets and possibly the pair- 
ing mechanism in the cuprate superconductors. 

First, we illustrated that the high resolution ARPES 
studies 3-4 from cuprates suggest the following. 

(a) The hole in the doped quantum antiferromagnet 
becomes a rather well-defined quasiparticlc and it is al- 
ways attached to a string of overturned spins which fol- 
lows the hole in its motion via processes of coherent pair 
spin-flip which lead to a reassignment of the hole birth 
site. 

(b) The high energy features and dispersion discussed 
in this paper, which are seen in ARPES, can be under- 
stood as "internal" excitations of the "spin-polaron" . In 
this paper, it is shown that these excitations are responsi- 
ble for the vanishing of the quasiparticle from the lowest 
string state at the T point and for the transfer of spectral 
weight to higher energy string excitations. In addition, 
they may be responsible for the high energy anomalous 
dispersion found in the most recent ARPES studies". 

Second, a simple model is postulated where only states 
of holes connected with unbroken strings of overturned 
spins are retained. This model gives a very good quan- 
titative account of both the quasiparticle dispersion and 
the higher energy string excitations. 

Third, motivated by both the fact that our simple 
model of string excitations can explain the results of the 



t — J model and the fact that the new features of ARPES 
studies can be interpreted as arising from such string ex- 
citations, we further consider the role of such string ex- 
citations in the pairing mechanism. We show that pairs 
of holes experience attractive pairing interaction due to 
exchange of such string excitations. Namely, when one 
of the holes creates such a string of overturned spins in 
order to be able to move, another hole takes advantage of 
the created path of overturned spins and follows it, thus, 
restoring part of the destroyed antifcrromagnetic order. 

The proposed pairing mechanism could be superior to 
those of models which are based on exchange of antifer- 
romagnetic spin fluctuations for the following two main 
reasons: (i) As we discussed in the previous section, our 
model does not require long range antiferromagnetic or- 
der but rather short-range antiferromagnetic spin cor- 
relations, (ii) The pairing energy scale is much larger 
than the typical energy scale that a spin-fluctuation ex- 
change model would give, since the effective attraction in 
our string-exchange model arises from exchange of short- 
range high energy excitations. 
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